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The CO oxidation behaviors on single Au atom embedded in two-dimensional h-BN monolayer are 
investigated on the basis of first-principles calculations, quantum Born-Oppenheim molecular dynamic 
simulations (BOMD) and micro-kinetic analysis. We show that CO oxidation on h-BN monolayer support 
single gold atom prefers an unreported tri-molecular Eley-Rideal (E-R) reaction, where O2 molecule is 
activated by two pre-adsorbed CO molecules. The formed OCOAuOCO intermediate dissociates into two 
CO2 molecules synchronously, which is the rate-limiting step with an energy barrier of 0.47 eV. By using the 
micro-kinetic analysis, the CO oxidation following the tri-molecular E-R reaction pathway entails much 
higher reaction rate (1.43 X 10^ s"^) than that of bimolecular Langmuir-Hinshelwood (L-H) pathway 
(4.29 s"^). Further, the quantum BOMD simulation at the temperature of 300 K demonstrates the complete 
reaction process in real time. 



Single atomic metal catalyst anchored to support has become attractive partly for the maximum usage of 
noble metal atom and great potential to achieve high activity and selectivity \ In experiment, the high 
surface free energy of metal atom is one key issue, which usually induces an aggregation of metal atoms into 
big cluster on support. In the past decades, appropriate metal oxide supports, which have strong interaction with 
metal atom, have been selected and several techniques have been developed to achieve supported single metal 
atom for catalysis^ ^ Meanwhile, the decreased freedom of catalyst in supported single metal atom materials 
makes it possible for theoretical understanding their catalytic behaviours. 

In recent years, two-dimensional monolayer materials, such as graphene, have been used as promising catalyst 
supports for their large specific surface areas, novel electronic and thermal properties, and close contact with 
catalysts, which may play an important role in activity enhancement of catalysts^"^°. To achieve atomically 
dispersed metal on graphene, various vacancy defects have been introduced into carbon network of graphene, 
strengthening the interaction between the metal atom and graphene. Both theoretical and experimental works 
have indicated that the single metal atom supported by defective graphene present enhanced catalysis, for 
instance the CO oxidation and oxygen reduction reaction on Au, Cu, Fe, Pt, and Nb supported by graphene^ 
However, the thermal stability and chemical reactivity issues associated graphene may hinder its application. 
Alternatively, the h-BN monolayer is a structural analogue to graphene, but exhibits quite different properties^^"^\ 
It is well known that h-BN monolayer is a wide band gap semiconductor with high chemical and thermal stability. 
Moreover, boron vacancy defect can be prepared by electron bombardment method^^'^^ and single metal atom can 
be deposited by mass-selected soft-landing techniques, improved wet chemistry method, or atomic layer depos- 
ition method, as that on grahene and other substrates^'^'^^'^^. The high ionicity of B-N bonds and possible vacancy 
defects in h-BN monolayer present great opportunity to modify the catalytic activity of supported single metal 
atom catalyst, even superior to graphene. Recently, Fe and Au atom supported by h-BN monolayer have been 
predicted theoretically as a single metal catalysis for CO oxidation^^"^^. 
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Figure 1 | Optimized structures of Au adsorbed and embedded h-BN monolayer. The optimized structures of Au atom on a) perfect and defective h-BN 
monolayer with either b) B vacancy or c) N vacancy are displayed from both the top and side views, respectively, d) to f) display the lowest energy structure 
of CO, O2, and two CO molecules on Au/Vb-/z-BN. The pink, blue, red, grey and yellow spheres represent B, N, O, C and Au, respectively. 



CO oxidation on gold nanostructures at low temperature has 
been intensively investigated for both advancing fundamental 
chemistry and industrial applications, including bare and oxide 
supported single atom, clusters, one- and two-dimensional struc- 
tures^'^^'^^"^^. In these gold-based CO oxidation catalysis reaction, 
either gold nanostructure itself (in gas phase) or the perimeter 
interface between the gold nanostructure and oxide support pro- 
vides active sites for the O2 activation via bimolecular Langmuir- 
Hinshelwood (L-H) mechanism, Eley-Rideal (E-R) mechanism, or 
both of them^^. A tri-molecular L-H mechanism has also been pro- 
posed recently when triangular Au sites are available on surface of 
gold clusters, where the O2 activation is promoted by CO molecules 
adsorbed on the gold clusters^^. The suggested factors that are 
responsible for high catalytic activity includes quantum size effect, 
structural effect of low- coordinated Au sites, strain of gold, charging 
effect, or combined effect of the oxide support and gold 
nanostructures^'^^'^^"^ ^ . 

In this work, on the basis of first-principles calculations and 
quantum Born-Oppenheim molecular dynamic (BOMD) simu- 
lation, an unreported CO-promoted O2 activation and tri-molecular 
E-R CO oxidation are firstly discovered on atomic gold supported by 
h-BN monolayer with boron vacancy (B vacancy). O2 is activated by 
two adsorbed CO molecules, and the formed OCOAuOCO inter- 
mediate dissociates into two CO2 molecules synchronously with an 
energy barrier of 0.47 eV, which is the rate-limiting step. Compared 
with bimolecular L-H mechanism, the trimolecular E-R mechanism 
presents a reduced energy barrier for CO oxidation. In contrast, Au 
atom adsorbed on the N vacancy defect will be displaced by O2 
molecule. The estimated CO oxidation rate following the trimolecu- 
lar E-R mechanism exhibit a high value of 1.43 X 10^ s~\ The 
quantum BOMD simulation at the temperature of 300 K reproduces 
the entire reaction process. 

Results and discussion 

We firstly examined the adsorption of a single Au atom on the 
pristine and defective h-BN monolayer with either N or B vacancy 



defect, denoted as Au/Zz-BN, Au/V^-h-BN, and Au/Vb-Zz-BN. The 
optimized structures are plotted in Figure 1. For the perfect h-BN 
sheet, Au atom prefers to adsorb on the top site of either N or B atoms 
with the adsorption energies of 0.10 and 0.11 eV, respectively. The 
distance between Au and N atom is as large as 2.93 A, as shown in 
Figure la. Similar distance value can be found for the adsorption on 
the top site of B atom. The weak interaction between Au atom and 
perfect h-BN sheet suggests that the adsorbed Au atoms may be 
clustering or easily desorbed from the h-BN sheet. Thus, a single 
Au atom supported by perfect h-BN sheet is not practicable for 
catalysis applications at room temperature. In contrary, a single Au 
atom can be anchored on the vacancy defect of h-BN sheet with short 
Au-N or Au-B bond. As shown in Figure lb and Ic, the shortest Au- 
N and Au-B bond lengths are 2.07 and 2.25 A for boron and nitrogen 
vacancy defect, respectively. The calculated adsorption energies of 
Au on Vn-Z^-BN and V^-h-BN sheets indicate that Au atom can bind 
strongly with the adsorption energy of 3. 1 7 and 3.45 eV, respectively. 
These values are comparable to the average cohesive energy of Au 
bulk (3.78 eV per atom)^^, indicating that a well distribution of single 
Au atom in the vacancy defect of h-BN sheet is achievable in experi- 
ment. The Au/Vb-/z-BN is ferromagnetic with a total magnetic 
moment of 2 |iB per supercell. The spin charge density mainly dis- 
tribute on the nitrogen atoms at the vacancy defect and the adsorbed 
Au atom. On the contrary, the Au/V^-h-BN is nonmagnetic. By 
placing Au atom in the vacancy defect, about 0.10 and 0.46 charges 
are transferred from Au to the defective h-BN sheet with nitrogen 
and boron defect based on Hirshfeld charge analysis respectively. 

Although Au binds to the defective h-BN monolayer strongly, it is 
still necessary to investigate the thermal stability of single Au atom on 
the defect in the reaction, i.e. the adsorption competition between the 
gas molecules (CO and O2 molecule) and Au atom on the vacancy 
defect of h-BN sheet. In an attempt to answer this question, we 
examined the adsorption energy of CO and O2 molecule on the 
defective h-BN monolayer. Various adsorption configurations have 
been considered (See Supporting Information Figure SI), where the 
largest adsorption energies are 8.81 and 1.75 eV for O2 on the V^-h- 
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Figure 2 | Reaction pathways of CO on Au-embedded h-BN monolayer, (a) CO oxidation follows a bimolecular L-H mechanism, (b) CO oxidation 
follows a tri-molecular E-R mechanism All energies are given with respect to coadsorbed CO molecules on Au/VB-h-BN. 



BN and Vb-/z-BN, respectively. For CO molecule, the values are 3.97 
and 5.86 eV for Vn-Z^-BN and Y^-h-BN, respectively. The large 
adsorption energy difference between the O2 molecule (8.81 eV) 
and Au atom (3.17 eV) on Vn-Z^-BN imphes that the adsorbed Au 
atom may be easily displaced when O2 molecule presents. The first- 
principle based BOMD simulation at 300 K confirms this view 
(Movies SI and S2 in Supporting Information). During 5 ps simu- 
lation times, the O2 molecule in gas phase displaces the adsorbed Au 
atom from Vn-Z^-BN and occupies the defect, whereas Au/Vb-/z-BN 
system is stable for CO and O2 in the simulation time scale. 
Interestingly, the B vacancy defects in /z-BN monolayer are found 
to be preferably formed in experiment by electron bombardment and 
theoretical calculation predicts that in N-rich environment B 
vacancy is more stable than N vacancy^^, suggesting the single atomic 
gold supported by /z-BN sheet with boron vacancy defect is practic- 
able for CO oxidation. 

In the following section, we only examined the CO oxidation 
process on single atomic Au supported by /z-BN monolayer contain- 
ing B -vacancy defect. Figure Id and le display the most stable struc- 
tures of CO and O2 molecule on Au atom. Different adsorption 



configurations and corresponding adsorption energies are summar- 
ized in Figure S2 (See Supporting Information). CO molecule is 
adsorbed on the Au atom with a tilted C-Au bond, and O2 molecule 
forms peroxide structure with elongated O-O distance of 1.33 A. The 
adsorption energy of CO molecule (1.28 eV) on Au/Vb-/z-BN is dis- 
tinctly larger than that of O2 (0.78 eV), implying that the adsorption 
of CO is energetically preferred. When CO presents on Au, the 
adsorption energy of O2 on Au/Vb-/z-BN is only 0.08 eV, whereas 
the adsorption energy of the second CO molecule is 0.22 eV without 
including van der Waals (vdW) interaction correction. Test calcula- 
tions with vdW correction also confirm this conclusion. Considering 
the vdW correction, these values are 0.40 and 0.18 eV, respectively. 
Thus, two CO molecules prefer to be adsorbed on a single atomic Au 
supported by defective /z-BN monolayer with an average energy of 
0.75 eV per CO molecule, forming a V-shape O-C-Au-C-O structure, 
as shown in Figure If The adsorption of two CO molecules on sup- 
ported Au largely results from the charge states of gold, where the gold 
atom stiU possesses +0.41 e charge with one CO adsorbed. Note that 
the adsorption of CO molecules on Au doesn't weaken the interaction 
between Au atom and the /z-BN monolayer with boron vacancy 
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Figure 3 | BOMD simulation of CO oxidation on Au-embedded A-BN 
monolayer. The snapshots at 0.0, 0.9, 2.4 and 4.0 ps of BOMD simulation 
of CO oxidation on Au/Vb-/z-BN are displayed. The pink, blue, red, grey 
and yellow spheres represent B, N, O, C and Au, respectively. 

defect. The calculated adsorption energy of AuCO and Au(CO)2 spe- 
cies with Vg-Zz-BN monolayer are 4.00 and 3.22 eV, respectively, 
which are similar with that of Au atom on Vb-/z-BN monolayer 
(3.45 eV). Thus, the single atomic catalyst of Au/Vb-/z-BN is ther- 
mally stable even during the CO oxidation process. The calculated 
energy barrier of CO molecule kicking off Au(CO)2 specie from Yb-H- 
BN is about 1.37 eV, which confirming the stability of the catalyst. 

Further, we explore the possible reaction pathways of CO oxida- 
tion on the single atomic Au catalyst anchored to the B vacancy of h- 
BN sheet. Very recently, Gao et al. investigated the O2 activation and 
CO oxidation on h-BN surface supported gold atom and dimer by 
using density functional theory^^ They proposed that the h-BN sur- 
face is not inert to gold nanoparticles and the CO oxidation can 
proceed via either a two-step pathway, i.e. two CO2 molecules are 
formed independently, or a self-promotion pathway on a nitrogen 
vacancy (N vacancy), following the bimolecular L-H mechanism^^. 
At first, the traditional bimolecular L-H mechanism is examined in 
our calculations (PATH-BLH), as shown Figure 2a. The energy of 
initial structure (IS) is set as zero. After one CO molecule is adsorbed 
on Au/Vb-/z-BN, O2 molecule attaches to the Au atom to form 
OCOO intermediate (MS-I in PATH-BLH) with an energy barrier 
of 0.23 eV. The OCOO intermediate dissociates into CO2 and an 
adsorbed O atom with an energy barrier of 0.72 eV (TS2-I). Further, 
the adsorbed O atom and gas CO molecule forms a CO2 molecule 
with an energy barrier of 0.15 eV (TS3-I). However, as the adsorp- 
tion energy of second CO on the CO pre-adsorbed Au/Vb-/z-BN 
(0.22 eV) is larger than that of O2 molecule (0.08 eV), the population 
of CO-O2 co-adsorbed configuration (fco-02) is almost two orders 
less than that of CO-CO co-adsorbed configuration (fco-co)j which 
can be estimated as fco-02/fco-co ^ exp(-AE/RT) with AE = 0.14 eV 
at room temperature (T = 293 K). Thus, the process via PATH-BLH 
may play a minor role in CO oxidation. 

In the following, the possible reaction among two adsorbed CO 
molecules and gas O2 molecule is investigated. Note that the gas 
phase O2 molecule travels in a spin-triplet configuration up to the 
distance close to the surface where hybridization becomes signifi- 
cant^^ and the adiabatic theory where its spin state continuously 



changes as it approaches the surface can result in an artificially small 
barrier for O2 dissociation. To overcome this problem, in our calcu- 
lations, the whole system is fixed at triplet state to provide reliable 
energy barriers for the O2 activation process, which has been adopted 
in similar simulations previously^^. In other reaction process, the spin 
states of whole system are fully relaxed to obtain the minimum 
energy reaction pathway. 

Figure 2b displays the reaction pathway of gas O2 molecule to two 
adsorbed CO molecules on Au/Vb-/z-BN. Different from PATH- 
BLH, as two CO molecules co- adsorb on Au/Vb-/z-BN, O2 activation 
prefers an unreported tri-molecular E-R mechanism (PATH-TER), 
i.e. O2 molecule is activated by CO molecules directly. Two paths 
with comparable energy barriers are found for the O2 molecule 
approaching CO molecules as shown in Figure 2b. In the transition 
state of the first pathway (TSl-II), O2 molecule approaches one CO 
molecule with the nearest C-O distance of about 1.65 A and O-O 
bond length of 1.31 A. The energy barrier is about 0.15 eV. In TS2-II 
structure, however, O2 molecule approaches two C atoms simulta- 
neously with a relative low energy barrier of about 0.09 eV. The 
distance between C atoms and O2 molecule ranges from 2.68 to 
2.79 A and O-O distance is about 1.24 A. The energy barriers in these 
two pathways are relatively lower than that in PATH-BLH (0.20 eV). 
All the structures of ISl-II, IS2-II, TSl-II, TS2-II, MS-II in the first 
step are in the triplet states, both reaction paths precede on triplet 
potential energy surface. Both reaction paths lead to the formation of 
OOCAuCOO* intermediate (—1.50 eV) with pentagonal ring struc- 
ture (MS), where O-O distance is elongated to 1.49 A. Next, 
OOCAuCOO* intermediate dissociates to form two CO2 molecules 
via TS2. The O-O distance is 1 .90 A and the energy barrier is 0.47 eV. 
The dissociation of OOCAuCOO intermediate is the rate-limiting 
step in whole reaction. The frequency calculations are performed to 
confirm the searched transition structures (See Supporting 
Information Figure S3). Also, the information about charge transfer 
and magnetic moment of whole system are summarized in Table SI 
(See Supporting Information). 

To visualize the complete reaction in real-time, we performed a 
first-principles based BOMD simulations at the temperature of 300 K 
to investigate the CO oxidation on Au/Vb-/z-BN (Move S2 in 
Supporting Information). The total simulation time is 5 ps, where 
the entire reaction process is reproduced. Figure 3 displays the snap- 
shots of structure at 0.0, 0.9, 2.4, and 4.0 ps. Initially, the distance 
between O2 and adsorbed CO molecule is longer than 5.0 A. With the 
time elapses, O2 molecule approaches to the adsorbed CO molecules 
and OOCAuCOO intermediate forms at about 0.9 ps. At about 
2.4 ps, O-O bond breaks up and two CO2 molecules are observed 
at about 4.0 ps, where one CO2 molecule is in gas phase and the other 
one weakly binds to Au atom. A structural optimization on this 
structure reveals that two CO2 molecules escape spontaneously and 
Au/VB-h-BN restores its original structure. Moreover, a BOMD 
simulation with four CO and two O2 molecules in gas phase also 
reproduces this process (See Movie S3 in Supporting Information). 

On the basis of micro -kinetics method (See details in Supporting 
Information), the maximum reaction rates for PATH-BLH and 
PATH-TER are calculated to investigate the favorability of two reac- 
tion pathways. CO oxidation following the tri-molecular E-R mech- 
anism entails much higher reaction rate (1.43 X 10^ s"^) than that of 
bimolecular L-H mechanism (4.29 s~^), which again confirms that 
the tri-molecular E-R mechanism is more favourable for CO oxida- 
tion in the Au/Vb-/z-BN system. In PATH-BLH, the reaction rate is 
mainly limited by the extreme small adsorption energy of O2 mole- 
cule and the high reaction barrier (0.72 eV) for TS2-II in Figure 2a. 
However, in the PATH-TER, the promoting effect of co-adsorbed 
CO molecules enables the formation of CO2 with much lower reac- 
tion barrier (0.47 eV). 

Alternatively, the third pathway (PATH-III) may be possible 
based on PATH-BLH, demonstrated in Figure S3 (See Supporting 
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Figure 4 | Electronic structure analysis of CO oxidation on Au-embedded A-BN monolayer, a) The profile distributions of O2S LUMO and HOS of CO 
coadsorbed Au/Vb-/z-BN. The isosurface value is 0.25. b) to d) are the PDOS profiles for CO coadsorbed Au/Vb-/z-BN, TS2-II, and MS structures, 
respectively. The red and green dotted lines denoted the DOS of free O2 or CO molecules, respectively. The Fermi energy level is plotted with black dotted 
line. 



Information). Instead of dissociation into CO2 molecule and O atom 
with a high energy barrier (0.72 eV), the formed OCOO (MS-I in 
Figure 2b) may react with the second CO molecule to form 
OOCAuCOO intermediate (MS-II) again with a relative low energy 
barrier of 0.23 eV. The rate-limit step in PATH-III is still the dis- 
sociation of OOCAuCOO into two CO2 molecules, as that in PATH- 
TER, whereas more reaction steps and higher energy barrier for O2 
activation are required in PATH-III than those in PATH-TER. Note 
that the CO oxidation via neither PATH-BLH nor PATH-III has not 
been observed in our BOMD simulations. 

The CO promoted O2 activation and synchronous generation of 
two CO2 molecules present a new tri-molecular E-R reaction for CO 
oxidation on supported Au atom, which has not been reported prev- 
iously. This kind of CO promoted O2 activation on Au/Vb-/z-BN can 
be understood with the frontier molecular orbitals involved in the 
reaction between O2 and co-adsorbed CO molecules. It has been well 
known that the O2 activation is accomplished by injecting extra 
electrons in the lowest unoccupied anti-bonding state (n'^2p)y i-e. 
the lowest unoccupied molecular orbital (LUMO) of O2 molecule. 
As shown in Figure 4a, the symmetric match between the LUMO of 
O2 molecule and the highest occupies states (HOS) of CO co- 
adsorbed on Au/Vb-/z-BN promotes the O2 activation. To further 
explain the interaction between CO, O2 and Au/Vb-/z-BN, we plotted 
the partial density of states (PDOS) of CO co-adsorbed Au/Vb-/z-BN, 
the TS2-II and MS-II structures in Figure 4b, 4c, and 4d. The DOS of 



free CO and O2 molecules are plotted in dotted line for comparison. 
As shown in Figure 4b, the CO adsorption on Au atom mainly 
originates from the interaction between Au's 3d orbital and CO's 
5a and 27i;* orbitals. At TS2-II structure (Figure 4c), a small amount 
of charge is transferred to 02's 7i;*2p orbital, which can be identified 
from 02's PDOS contribution at the Fermi energy level in Figure 4c. 
It is clear that 02's 7i;*2p state is fully occupied for the major spin, 
while that for the minor spin is splitted and partially occupied. In the 
MS-II state, strongly hybridized states between co-adsorbed CO and 
O2 molecules can be found just below the Fermi energy level, as 
shown in Figure 4d. The highest occupied states (HOS) and the 
lowest unoccupied states (LUS) of MS-II state are shown in Figure S5. 

At last, the CO oxidation on Au dimer supported by /z-BN mono- 
layer with B vacancy, i.e. Au2/Vb-/z-BN, was also examined. The 
adsorption structures and energies are summarized in Figure S6 
(See Supporting Information). Au dimer is trapped in the vacancy 
with a tilted Au-Au bond to the /z-BN monolayer surface and the 
average adsorption energy is 2.53 eV per Au atom. This value is 
remarkably smaller than that of a single Au atom, indicating a well 
distribution of Au atom on the defective /z-BN monolayer. Up to 
three CO molecules can be adsorbed on Au2 with one on the Au atom 
at the vacancy (Au-1) and two molecules on the other Au atom (Au- 
2). The calculated adsorption energies for three CO molecules are 
1.30, 1.11, and 0.25 eV in sequence. Figure 5 illustrates the reaction 
pathway for CO oxidation. The O2 molecule reacts with the CO 
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Figure 5 | CO oxidation on Au2-embedded A-BN monolayer. Reaction 
pathway for CO oxidation on Au2/VB-h-BN. The energy given in figures is 
the energy relative to initial state. 

molecule on Au-1 atom, forming an intermediate OCOO structure 
with the energy barrier of 0.17 eV at the first step. Next, the OCOO 
intermediate reacts with one CO molecule on Au-2 and forms 
OCOOCO ring structure with an energy barrier of 0.10 eV. At the 
last step, the OCOOCO intermediate dissociates into two CO2 mole- 
cules with an energy barrier of 0.37 eV, which is 0.10 eV lower than 
that on Au/Vb-/z-BN. 

In conclusion, the CO oxidation on single Au atom anchored to 
the h-BN monolayer with boron vacancy defect was examined by 
using first-principles calculations and quantum BOMD simulation. 
We found that Au atom on N vacancy will be displaced by gas CO 
and O2 molecule, while much stable on B vacancy. The CO oxidation 
on Au atom supported by /z-BN monolayer with B vacancy prefers an 
unreported tri-molecular E-R mechanism, where O2 is activated by 
two adsorbed CO molecules with an energy barrier of 0.09 eV. The 
formed OOCAuCOO intermediate dissociates into two CO2 mole- 
cules synchronously, which is the rate-limiting step with an energy 
barrier of 0.47 eV. Compared with traditional bimolecular L-H 
mechanism, trimolecular E-R mechanism largely demands lower 
reaction energy barrier. By using the micro -kinetic analysis, the 
CO oxidation following the tri-molecular E-R mechanism entails 
much higher reaction rate (1.43 X 10^ s"^) than that of bimolecular 
L-H mechanism (4.29 s~^). The quantum BOMD simulations repro- 
duce the whole reaction. Our results demonstrated a potential for CO 
oxidation on h-BN monolayer supported atomic gold. 

Methods 

All calculations are performed by using the linear combination of atomic orbital and 
spin-polarized density functional theory method, implemented in DMoP package 
(DM0I3 5.5 is a density functional theory quantum mechanical package available 
from Accelrys Software Inc.)^''"*^° The generalized gradient approximation with 
Perdew-Burke-Ernzerhof (PBE) functional form with effective core potential, and 
double numerical basis set with polarization function are adopted in the calcula- 
tions'^^'^^l The real-space global cutoff radius is set to be 4.5 A. A tetragonal supercell 
containing 6X6 unit cells of /z-BN monolayer is used and the thickness of vacuum 
layer is chosen as 20.52 A. Only F-kpoint is used to describe the Brillouin zone for 
geometric optimization and 6X6X1 kpoint for electronic structures calculations^"^. 
The smearing parameter is set to be 0.001 Hartree. Test calculation indicates that the 
results don't change when the smearing parameter is set zero. The structures of TS are 
determined with synchronous transit methods and confirmed with frequency cal- 
culations. The constant temperature and constant volume quantum BOMD simu- 
lation is performed by using the Nose-Hoover method, in which the kinetic energy 
fluctuation of the thermostat variable is controlled by coupling it to another ther- 
mostat variable and the temperature of the system is maintained at 300 K'^^ '^^. The 
time step is 1.0 fs and the total simulation time is 5.0 ps. Note that the spin crossover 
effect is not considered in the calculation, however, the spin stats of reactants, 
transition states, intermediates, and products are fully examined to obtain the lowest 



energy structures. The CO oxidation rates are estimated by using the micro-kinetic 
analysis (See details in Method in Supporting Information). 
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